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O . Abstract 

On 

On ' We have formulated a self-consistent model of freeze-out on an arbitrary hyper- 

surface. It conserves energy and momentum across the discontinuity between ideal 
fluid and the gas of free particles. Energy and momentum of those free particles 

Q ' have non-equilibrium values that could be a signal for the formation of hot and 

dense matter in heavy ion collisions. 
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1. Introduction 

In this paper we present our model of freeze-out for relativistic hydrodynamics. It 
solves the problem where and how to stop solving the hydrodynamical equations in such a 
way that it is consistent with the conservation of energy and momentum. When one solves 
the hydrodynamical equations one has to terminate the solution on the boundary with 
the vacuum. Also, one has to stop solving hydrodynamics at some freeze-out temperature 
or baryonic density. However, in most cases one continues to solve the equations for 
such values of these quantities, for which hydrodynamics is no longer valid (see p], H). 
There are several papers || || on a shock-like treatment of the freeze-out problem. 
However, in our opinion these approaches are ad hoc, since the existence of such a shock 
is postulated and not obtained as a result of the equations of motion. 

It should be noted that this is only part of the trouble. Another one comes from 
the problem of calculating particle spectra on a time-like hypersurface. For a space-like 
hypersurface of freeze-out the correct answer for the spectra of particles is given by the 
formula of Cooper and Frye ||. However, one cannot use that formula for time-like 
hypersurfaces, since it leads to negative numbers of particles. This is due to the fact that 
it was obtained only for the space-like case, where the decay of one element of gas does 
not affect the decay of adjacent elements. 

Here we present our model which is based on the conservation laws of energy and 
momentum between fluid and gas of free particles. The paper is organized as follows: 
the second part contains a short derivation of the momentum spectra for the gas of free 
particles on a time-like hypersurface, which was also obtained recently in H; the third 
part contains the derivation of the conservation laws; some useful formulae are derived 
there; and in the final part we discuss possible effects on observables. 

We also hope that the model suggested here is not only of academic interest, but 
will be used by other researchers to solve the equations of relativistic hydrodynamics in 
practical calculations. 



2. Decay of the gas of free particles 

In order to obtain the particle spectra for the gas, we will use the method derived 
by Gorenstein and Sinyukov ||. In this paper we shall deduce the equations for the 
case of one-dimensional hydrodynamical motion, but the final result will not depend on 
this assumption, since it will be written in covariant form. Suppose there is a boundary 
between fluid and gas. Let us consider the decay of a small element Ax of the gas of 
free particles in its rest frame. The gas is supposed to be located in the left hemisphere 



and to have the freeze-out temperature T = T*. (We suppose that the derivative to the 
freeze-out hypersurface v ff in the t — x plane is positive). We note that this frame is the 
rest frame of the gas before decay. Hereafter we shall call it the reference frame of the 
gas. Suppose the particles in the gas have an equilibrium distribution function <f> (§£)• 

First we consider the contribution from particles with negative momenta that leave 
the element Ax (see Fig.l) 
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where p± is the transverse momentum of the particle, and AS± the transverse size of the 
element. The second contribution is given by particles with negative momenta from the 
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Finally, the third contribution comes from particles with positive momenta from the 
element Ax — %-At. However, those particles will cross the freeze-out hypersurface only 
if their velocity is smaller than the derivative to the hypersurface v a in the t — x plane. 
Thus, the third term reads as follows (see Fig. 2): 
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After some simple algebra one obtains the formula for the spectrum of the gas of free 
particles 
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As will be shown below, the modification of the spectrum due to the G-function will lead 
to an energy-momentum tensor that differs from the equilibrium case. 

Our last step is to write the formula for the spectrum in a fully relativistic form. For 
that we need to change the energy of the particles appearing in the distribution function 
in the reference frame of the gas to the product of the four-vectors of momentum and 
hydrodynamical velocity, p^u^ 1 , and change the integration over the hypersurface of freeze- 
out to the product of the four- vectors of momentum and normal vector to the freeze-out 
hypersurface, p^da^ 1 . Finally, we have 



po J^t = K^-) p " d(xueippdap) ' (5) 

where the vector da^ = (v a , —l)dt is the normal vector to the freeze-out hypersurface in 
the left hemisphere. It is, however, easy to check that the above formula is valid for the 
right hemisphere as well. 

The meaning of this nice result is that it is the formula of Cooper and Frye || , but 
without negative particle numbers! 

It is easy to see that for a space-like hypersurface, where v a > 1, the above expression 
gives the result obtained by Cooper and Frye 0. 

The energy-momentum tensor of free particles in the reference frame of the gas has 
the following form: 

It is easy to calculate this tensor for the case of noninteracting massless particles; it has 
the form: 

r 2 °v) = e (n 1 ^, (7) 

T 2 > CT ) = e (T*) ^^ , (8) 

r 2 >,) = e(T*)^±I, (9) 

where e is the usual energy density. The above result is valid for the left hemisphere. The 
corresponding formulae for the right hemisphere can be obtained in the same way. Now 
one can see that for the case v a = 1 the above expressions give the usual formula for the 
ideal gas. We hope that this unusual behavior of energy and momentum can be found in 
experimental pion spectra. 

In the paper by Sinyukov || another model for freeze-out was suggested. He considered 
the decay from a box and did not account for the additional contributions from the 
intrinsic volume of the gas, namely from the element —2-At. Due to that, the energy- 
momentum tensor obtained in his paper is not symmetric! It means that the orbital 
momentum of the considered system is not conserved! 

On the other hand, he considered an ad hoc model for deflagration from very hot pionic 
matter into gas of free particles. We do not think that such a simple picture corresponds 



to the real situation in heavy ion collisions, since the solution of the hydrodynamical 
equations fl], || |9| does not exhibit shock-like transitions in the expansion of hot and 
dense pionic matter. 

Now we would like to clarify an important question: "What is the difference between 
an ideal fluid and the gas of free particles at freeze-out temperature?" It seems that the 
main difference is that they have different values of the cross-section. Due to that, there 
are collisions in the fluid which lead to thermodynamical equilibrium. In contrast, there 
are about no collisions in the gas of free particles (and we shall neglect them completely), 
because the cross-section for collisions is very small. Of course, fluid and gas have some- 
what different values of temperature, but due to the fact that the mean free-path strongly 
depends on the temperature, this difference should be small. It was found [|l(| that for 
pions the mean free-path depends on the fifth power of the inverse temperature: A ~ ^^-. 
Thus, the difference between the temperatures of the fluid and the gas should be small, 
but due to the strong dependence of the collision cross-section on the temperature, their 
mean free-paths should be very different. 

Now we would like to formulate a more physical concept for freeze-out, based on the 
conservation laws on the discontinuity between ideal fluid and the gas of free particles. 



3. Conservation laws on the surface between fluid and gas 

We would like to consider a simple model of freeze-out assuming that the system consists 
of the fluid and the gas of free particles. We do not take into account the transition region 
between them. Then, the total energy-momentum tensor of the system is as follows 

t^ = rre(r - t*) + t£ u q(t; - t) , (10) 

where index 1 corresponds to the fluid, and index 2 to the gas of free particles, T is the 
temperature, Tj* and T 2 * are the freeze-out temperatures for fluid and gas, respectively. 
We assume that those temperatures satisfy T* > T 2 *. 
The equations of motion have the form: 

dpT 1 ™ = . (11) 

Taking derivatives and using the equations for the evolution of fluid and gas, 

a M rr = 0, a = 1,2, (12) 

one finds that the terms with the delta-functions must vanish: 



T^d^T* = T^d^T* . (13) 

Thus, we have obtained two equations (for the 1+1- dimensional case, but the case of 
3+1 dimensions gives the same result). However, we have to add to this equation the 
condition that the derivatives of the temperatures of fluid and gas are equal, since there 
is only one hyper surf ace of freeze-out: 

d^T* = d^r; . (14) 

Using this condition and dividing the first equation above by the second one and 
canceling derivatives of T*, we get a new equation: 

( T ii _ t] 1 ) (r°° - T 2 00 ) = (T™ - T 2 10 ) (T* 1 - T 2 01 ) , (15) 

or 

Det(TF -Tf) = . (16) 

However, there is one problem. It seems that we know the energy-momentum tensor 
of the gas only in its reference frame, but the energy-momentum tensor of the fluid is 
known in the laboratory system. On the other hand, we do not know the velocity of the 
gas of free particles in the laboratory system. Fortunately, this problem can be easily 
solved. 

Let us suppose that the fluid with temperature T = T x * has a velocity v\ in the 
laboratory system, and the gas has a velocity t>2 in the same system. Then, the velocity 
of the fluid in the reference frame of the gas is 

Vrel = ~ • (17) 

1 - V^ 2 

Now the energy- momentum tensors of the fluid Ti V (v re i, T*) and the gas T^^iv, T 2 *) are 
known in the reference frame of the gas. And we have one equation that was obtained 
before: 

Det (Tf>rez, T*) - 7T K, T*)) = . (18) 

Fortunately, we have one more equation. Due to the definition of the freeze-out hyper- 
surface: 



T(x,i) = T 2 * (19) 
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and the fact that the velocity v ff is the derivative of the implicit function T(x, t) = T 2 
one gets the following equation 

/ cfe\ _doTl n°(v reh T*)-T™(v a ,TZ) 

Va \dt) T=T , d 1 T 2 * T™(v rel ,T?)-T™(v ff ,T 2 *)- [ Uj 

Another way to obtain the above result is as follows. One has to remember that the 
derivatives of the freeze-out temperature T 2 are the components of the normal line to the 
hypersurface T(x,t) = T 2 . Then, energy-momentum conservation reads 



Tl\v reU Tl)-Tl\v reh T*)v 
T^(v rel ,T*)-T^(v rel ,T*)v c 

Thus, we have two equations for three unknowns Tj*, v re i and v ff , since we suppose that 
the freeze-out temperature for the gas of free particles is known. 

From the above equations one finds the freeze-out hypersurface in the reference frame 
of the gas. However, this hypersurface is still unknown in the laboratory frame, because 
the above equations do not fix the velocity of the gas in the laboratory frame. In order 
to find the freeze-out hypersurface in the laboratory frame, one has to solve the hydrody- 
namical equations. Thus, the problem is solved. 

We have to add the following. The equations 

TfdpTZ = TfdpTZ (23) 

are the conservation laws of energy and momentum on the discontinuity between the fluid 
and the gas of free particles, since the derivatives of the freeze-out temperature T* give the 
normal vector to the hypersurface. However, one cannot obtain them from the solution 
of the equations of hydrodynamics for the fluid. Suppose one knows the freeze-out point 
T(x, t) = T* at a fixed time. In order to find this point at the next moment of time, one 
has to solve the hydrodynamical equations for unphysical temperatures, namely, T < T*\ 
Thus, from the above equations one has to find the velocities v ff and v re i in the reference 
frame of the gas. Since the velocity of the fluid in the laboratory system v\ is known from 
the solution of the hydrodynamical equations, one makes a Lorenz transformation and 
finds the velocity of the gas in the laboratory system: 

vi-v rel 

v 2 = . 24 

1 - ViV rei 



Using this velocity one can find the derivative of the freeze-out hypersurface in the labo- 
ratory frame by a Lorenz transformation. 

In the case of baryon-rich matter, the freeze-out condition is as follows 

n(x,t)=n*, (25) 

and the thermodynamical functions depend on temperature and chemical potential. In 
this case we have one more hydrodynamical equation which is the conservation law of 
baryonic charge, and one more equation for the baryonic densities between fluid and the 
gas of free particles. 



4. Results and discussion 

We have developed a freeze-out model for relativistic hydrodynamics. We have ob- 
tained the conservation laws on the discontinuity between ideal fluid and the gas of free 
particles. We have derived an expression for the energy-momentum tensor of the gas of 
free particles and its momentum distribution function. It is important to emphasize that 
this tensor differs from the equilibrium one. Thus, one can hope to find this unusual be- 
havior in experiments, and this can give detailed information about the freeze-out process 
in heavy ion collisions. 

One possible effect is related to measuring the ratio of the total (longitudinal) momen- 
tum density in the backward hemisphere to the zero component of the baryonic charge 
flux. In the reference frame of the gas, that component reads (for the left hemisphere) 

iV° = n(:r»i±^. (26) 

If it is possible to neglect the contribution of the space-like part of the freeze-out hyper- 
surface, then that ratio will depend on the thermodynamical quantities and an average 
value of v a on the time-like part of the freeze-out hypersurface. Thus, measuring this 
ratio for the gas of free particles one can find the average value of v a . We hope that it is 
possible to do so for asymmetric heavy ion collisions. For example, the S + Au reaction 



that was studied in [|II]] indicates that the space-like part of the freeze-out hypersurface 
is small and that the velocities on it are approximately constant. Therefore, finding the 
suggested ratio can give important information about the freeze-out process of hot and 
dense hadronic matter. 
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Appendix A 

In this Appendix we would like to discuss the question: How to solve the hydrody- 
namical equations together with the above equations on the freeze-out hypersurface. In 
fact we shall only show that the suggested system of equations for the shock is consistent 
with the hydrodynamical equations, so that they can be solved simultaneously. 

A. Suppose our system consists of pions only. In this case thermodynamical quantities 
depend only on the temperature. Suppose we have two hydrodynamical equations for the 
fluid that can be written in the form 

F a {€ii,Pii,v u ,Xi) = 0, a = 1,2 , (27) 

with unknown values of energy density en, pressure pn, velocity vu, and at the coordinate 
point Xi. Now let us consider the evolution of the discontinuity between fluid and the gas of 
free particles. First, let us count the unknowns. They are 7\, v r ^ , v±, v 2 (instead of v re i) 
and xl hock . The last one we can change to the velocity of the shock in the laboratory 
system: v^ ab . Thus, we have 5 unknowns. On the other hand we have precisely 5 
equations. They are: two hydrodynamical equations for the fluid, two conservation laws 
on the discontinuity between fluid and gas and the relation between the velocities in 
different frames: 

1 + V a V 2 

Thus, solving this system of equations one finds the trajectory of the freeze-out hypersur- 
face, temperature of the fluid and velocities of fluid and gas. 

B. Now let us consider the case of baryonic matter. Then one has more unknowns and 
more equations. Let us count them. 

The unknown quantities are: three thermodynamical functions of the fluid - e\, p\, 
rii, two thermodynamical functions for the gas - T 2 , fj, 2 , the derivative to the freeze-out 
hypersurface in the reference frame of the gas - i£ e ' , the velocities of the fluid and gas 
in the laboratory system v\, v 2 , and the derivative to the freeze-out hypersurface in the 

9 



laboratory system v^ ab . Thus, one has 9 unknowns. At the same time there are also 9 
equations: 3 hydrodynamical equations for the fluid (two for energy and momentum, one 
for the baryonic charge), 3 equations on the discontinuity between fluid and gas (two for 
energy and momentum, one for the baryonic charge), the equation of state for the fluid 
that relates its thermodynamical functions e±, pi, rii, the relation between the velocity of 
the shock in the reference frame of the gas and in the laboratory frame, and the equation 
of freeze-out for the gas of free particles n 2 (x,t) = n*. Now one has to solve the above 
system of equations and find all unknowns. 
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Fig.l. Decay of the element A x = AD. 
Particles with negative momenta are 
leaving elements AD and DE = - p/p^ A t . 
The passage of particles is indicated by lines 



with arrows. 
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Fig.2. Decay of the element A x = AD. 
Particles with positive momenta are 
leaving elements AF = A x - p/p At . 
The passage of particles is indicated by lines 



with arrows. 
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